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ABSTRACT 

The determination of the arrival direction of gamma-induced air showers in the 
ARGO-YBJ experiment has been investigated using different algorithms in the 
reconstruction procedure. The calculation has been performed, as a function of 
pad multiplicity, for different primary energies. The performance of a conical 
correction to the shower front with a suitable fixed cone slope is discussed. 



1 Introduction 

The ARGO-YBJ detector, currently under construction at the Yangbajing Lab- 
oratory (P.R. China, 4300 m a.s.l.), is a full coverage array of dimensions 
^ 74 X 78 realized with a single layer of RPCs. The area surrounding 
this central detector (carpet), up to ~ 100 x 110 m^, is partially (~ 50%) 



instrumented with RPCs. The basic element is the logical pad (56 x 62 cm^) 
which defines the time and space granularity of the detector. 

In a search for cosmic point-like sources the main problem is the back- 
ground due to charged cosmic rays, so a good angular resolution (measurement 
of the arrival direction) is necessary. In order to determine the primary direc- 
tion, the shower front has to be reconstructed from measurement of the time 
and position of the fired detectors. This calculation is usually performed for 
internal events only. In fact, the direction of showers with the core outside 
the detector is in general badly reconstructed due to the cone-like shape of the 
shower front and the unknown core position. 

In Fig. ^ is shown an exaggerated shower with two scenarios: case A 
where the shower core strikes at the center of the ARGO-YBJ detector, and 
case B where the core hits outside the carpet. The primary energy in case 
B must be larger than that in case A to be able to trigger the array with 
the same detector response (number of fired pads). In case A, the direction 
obtained from the plane fitting is close to the true shower direction. On the 
other hand, in case B, the carpet sees a tilted portion of the same air shower, 
and in turn the fitted direction acquires a pointing error and deviates from the 
true shower direction. 

In this paper we present different algorithms to reconstruct the shower 
direction. The calculations have been performed for different primary energies. 
The effect of a lead converter sheet on the detector is also investigated. The 
performance of a conical correction to the shower front with a suitable fixed 
cone slope is discussed. 

2 Reconstruction procedure of the shower direction 

The arrival time of the earliest particles provides the information concerning the 
time profile of the shower front. From MC simulations (see Fig. 15 of ref.BJ) it 
appears that the shower front assumes a parabolic shape more pronounced for 
the electron component than for the photon component. Near the core charged 
particles form, at a first approximation, a plane disk. 

The thickness of an EAS is defined by the distribution of arrival times of 
the shower particles which are delayed relative to the prompt particles which 
make up the shower front. Therefore, the shower thickness is determined by 
the lower-energy particles in a shower. This means that the shower thickness 




Figure 1: Core inside or outside the 
carpet can make the difference. The 
true shower core hits under the central 
arrow. 



Figure 2: Distributions of temporal 
residuals calculated for 1 TeV verti- 
cal ^-induced showers after a different 
number of plane fit iterations. 



is the smallest near the core and increases with the increasing distance from 
the shower axis position © . Hence, even in the case of vertical incidence the 
particles arrivals are not simultaneous. 

In addition to the curvature, the combined effects of the finite size and 
time resolution of the detector make the shower apparently non-planar: since 
the shower plane has a thickness and the arrival time of the first particle is 
measured by the TDCs, the larger the number of particles, the earlier the 
measurement (on average). 

The usual method for the reconstruction of the shower direction is a 
fit to the recorded arrival times ti by minimization of 

X^ = ^^/-i.)' (1) 

i 

where the sum includes all pads with a time signal. Usually the function / 
describes a plane, a cone with a fixed cone slope, or a plane with curvature 
corrections depending on r and N . The weights w are generally chosen to be 
an empirical function of the number N of particles registered in a detector. 

An improvement to this scheme can be achieved by excluding from the 
analysis the time values belonging to the non-gaussian tails of the arrival time 



distributions with successive minimizations for each shower EJ. 

In this paper we will investigate the performance of the ^ fit reconstruct- 
ing the primary direction of the showers sampled by the ARGO-YBJ carpet by 
means of the following iterative procedure: 

a) Planar fit 

In the first step we recover the shower direction cosines {Zp, rrip} by means 
of an unweighted planar fit performed minimizing the function 

^{ns^) = \ ^{ZaJi + mj/j + nzi + c{U - ^o)}^ (2) 

i 

The sum includes all pads with a time signal f^, c is the velocity of light, 
(xi, Hi, Zi) are the coordinates of the central position of the i-th pad. The 
parameters of the fit are the time offset and the direction cosines m. 
After each minimization the time signals which deviate more than K ■ a 
from the fitting function are excluded from further analysis and the fit 
is iterated until all times do not verify this condition or the remaining 
hits number is < 5 (in this case the event is discarded). Here a is the 
standard deviation of the time distribution around the fitted plane (i.e., 



nZ^)- III principle, the value of K could depend on the features of 
the reconstructed showers as well as on the experimental condition (pad 
dimension, shaping of the signals, multihit capability, etc.). 

b) Conical correction 

With this first determination of the arrival direction {Zp, mp} we calculate 
the conical correction 

SU = -- R, (3) 
c 

where 

R, = ^ (Axf + Ay2) _ (Aa;,/p + Ay.mp)^ (4) 

and Axi ~ x-i ^ Xc, Ay^ — Hi — Vc are the pad distances from the shower 
core position {xc, yd- Then we correct the experimental time values ti 

ti~^t[^ ti - Sti{lp, nip) (5) 



In this approach the slope parameter a is not a fit parameter but it is 
fixed to a value that improves the angular resolution. 
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Figure 3: Opening angle ■(/'yo '^'■5. num- 
ber of iterations for different values 
of the K parameter. The upper plot 
refers to 1 TeV ^-induced showers, 
the lower one to 500 GeV (Nhit = 
20- 50/ 



Figure 4: Opening angle ipjo vs. num- 
ber of iterations for different values 
of the K parameter. The upper plot 
refers to 1 TeV ^-induced showers 
(Nhit 100 - 200), the lower one to 
5 TeV (Nhit = 500- 600;. 



c) New planar fit 

With the time values t'^ we reconstruct a new shower direction by means 
of a further planar fit. 

In Fig. Elthe reduction of the non-gaussian tailes as a function of the successive 
plane fit iterations is shown. After few iterations the time residuals distribution 
assume a symmetrical shape (dotted histograms). 

In order to optimize this procedure we have studied the resulting angular 
resolution as a function of the number of planar fit iterations and of the values 
for K. To estimate the performance of the algorithms we have used the V'TO 
parameter, defined as the value in the distribution of the angle between the 
true and the reconstructed directions within which 71.5 % of the events are 
contained. The angular resolution ag is given by ae = 1.58 • V^to- In Fig. Elthe 
behaviour of the opening angle ^70 as a function of the number of iterations 
is shown, for different choices of K. The upper plot refers to 1 TeV 7-induced 
showers, the lower to 500 GeV 7-induced showers. In these calculations, all 
showers have been sampled with the core at the center of the carpet and we 



have selected events with a pad multiphcity Nhu in the range 20 - 50. The 
same plots for higher multiplicities are shown is Fig. 0] No converter has been 
added on the RPCs. 

From these figures it results: 

• The angular resolution improves substantially with successive iterations. 
Indeed, the successive minimizations, as expected, reject the farthest tem- 
poral hits from the "theoretical plane" (see Fig. This behaviour is 
qualitatively independent of the primary energy, the pad multiplicity and 
the values of K. 

• After a few iterations the improvement stops and, for high multiplicities, 
the angular resolution worsens (see Fig. ^J. This feature is due to the 
large fraction of discarded temporal hits lying far from the shower core. 
As a consequence, we must fix a suitable maximum number of iterations. 

A systematic study "^^ performed for different energies and zenithal angles, 
sampling the shower cores over a wide area, led us to the following best tuning 
of the planar fit procedure: K = 1.5 — 2 with the maximum number of iterations 
in the range 4-6. 

A further improvement of the angular resolution can be obtained by fixing 
the cone slope of the successive conical correction to the value a = 0.03 ns/m. 
As an example, in Fig. [S] we show the parameter ?/'70 a function of pad 
multiplicity for different values of the cone slope a. The calculations refer to 
vertical 1 TeV 7-induced showers with the core at the center of the carpet. As 
can be seen, the value a — 0.03 ns/m improves the angular resolution for all 
the investigated multiplicities (see Fig. 

We also investigated the performance of a conical fit in which the cone 
slope is a free parameter of the minimization El . The angular resolution cal- 
culated with fixed a gives better results because the minimization tends to 
adapte the cone slope to the temporal hits, thus introducing a larger error in 
the direction recontruction for events with the core outside the carpet or in its 
boundaries EJ. 

2.1 The effect of a lead converter on the angular resolution 

The consequences of placing a thin sheet of converter above the detector are, 
qualitatively: (1) absorbtion of low energy electrons (and photons) which no 
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Figure 5: The parameter ip^Q as a 
function of pad multiplicity for differ- 
ent values of the conical parameter a 
in vertical 1 TeV ^-induced showers. 
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Figure 6: Opening angle ipjo as a 
function of the conical parameter a 
for vertical 1 TeV ^-induced showers 
with a selected pad multiplicity Nhit = 
200- 250. 



longer contribute to the time signal; (2) multiplication process of high-energy 
electrons and photons which produces a signal enhancement. The lead con- 
verter thus reduces the temporal fluctuations: the contributions gained are 
concentrated near the ideal time profile since the high energy particles travel 
near the shower front while those lost tend to lag far behind. 

In order to study the effect of the converter on the angular resolution we 
have simulated, via a GEANT3-based code, a 0.5 cm lead sheet located 5 cm 
above the RPCs. The results discussed in the previous section are qualitatively 
reproduced also when the converter is added for what concerns the dependence 
of the angular resolution on the number of iterations of the planar fit. As 
expected, the improvement of the showers temporal profile implies a fewer 
number of iterations in the iterative procedure: the best tuning of the planar fit 
is achieved after only 3 iterations, with K = 1.5 . In Fig. [7|we show the effect of 
a 0.5 cm lead sheet on the angular resolution by comparing the "070 parameter 
as a function of the fired pads with and without the lead. The calculations 
refer to 1 TeV 7-induced showers with a fixed zenithal angle 9 = 20° and with 
the core sampled on a 50 x 50 m^ area. In the upper plot we show the values 
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Figure 7: The parameter tJjiq as 
a function of pad multiplicity after 
the planar fit (upper plot) and after 
the conical correction (lower plot) for 
1 TeV ^-induced showers (6 — 2Q° ). 
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Figure 8: The parameter ipia o,s a 
function of pad multiplicity for 7- 
induced showers. We compare the 
opening angle calculated with different 
algorithms. 



calculated after the best tuned iterative plane fit: 5 iterations without lead and 
3 iterations with the converter. In the lower part of the figure "070 is calculated 
after the conical correction. The improvement is clearly evident. 

3 Angular resolution of the ARGO-YBJ detector 

In oder to study the angular resolution of the ARGO-YBJ detector we have 
simulated, via the Corsika codelS*, 7-induced showers with a Crab- like spec- 
trum (~ and proton events with ^ i?^^-^^, both ranging from 100 GeV 
to 50 TeV. The 7-rays have been simulated for different zenith angles (< 40°), 
following the daily path of the source in the sky. The detector response has 
been simulated via a GEANT3-based code. A 0.5 cm lead sheet have been con- 
sidered 5 cm above the RPCs. The core positions have been randomly sampled 
in an area, energy-dependent, large up to 800 x 800 m^ and centered on the 
detector. 

In Fig. |S1 we compare the dependence of the opening angle 7/170 on pad 
multiplicity for the different algorithms discussed above. The 7-induced show- 
ers are reconstructed inside a fiducial area similar to the carpet dimensions 
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Figure 9: The parameter ip^Q as a 
function of pad multiplicity for 7- 
induced showers after the conical cor- 
rection. 
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Figure 10: The parameter tpjQ as a 
function of pad multiplicity for 7 and 
proton showers. 



{A fid = 80 X 80 m^). The shower core positions are reconstructed by means 
of the Maximum Likehhood Method appHed on the RPCs (the so-called 'LLF2 
method' IS'). The conical correction needs a good shower core position resolu- 
tion, so for very low multiplicities the planar fit is more performant because of 
the high contamination of external events. However, when suitable "cuts" ^ 
are applied and a high fraction of external showers is rejected, the conical cor- 
rection improves the planar fit. In Fig. |^ we show the parameter & 
function of pad multiplicity for 7 showers reconstructed inside ^/id with and 
without the application of the cuts. For comparison, the '070 values for the 
events with the sampled core truly inside A fid are also plotted (circles). As 
can be seen, for low multiplicities {Nhn < 100), the rejection of external events 
is crucial to improve the angular resolution. 

In Fig. 1101 a comparison between the angular resolutions obtained after 
the cuts for gamma and proton showers is shown. Due to the larger lateral 
spread of particles in proton showers, the angular resolution for protons is 
worse than that for photons. 



4 Conclusions 



In this paper we described an algorithm optimized to reconstruct the direction 
of the showers detected by the ARGO-YBJ cxpcriniont. A conical correction 
with a cone slope parameter fixed to the value a. = 0.03 ns/m is proposed and 
discussed. The resulting opening angle i/'TO is better than 0.3° for 7-induced 
events which fire more than 600 pads on the carpet. 

For Nfiit > 100 (when the dependence on the fiducial area is reduced) the 
angular resolution of the ARGO-YBJ experiment, for 7-induccd events with 
a Crab-like energy spectra, can be described by the equation ip-jo = 0.046 + 
6.666/V^. 

References 

1. G. Di Sciascio, B. D'Ettorre PiazzoH, M. lacovacci, Astrop. Phys. 6 (1997) 
313. 

2. D. Campana, B. D'Ettorre Piazzoli and G. Di Sciascio, Nucl. Instr. and 

Mcth. A344 (1994) 250. 

H. Krawczynski et al., Nucl. Instr. and Meth. A383 (1996) 431. 

3. E. Rossi, University Degree Thesis, Naples (2002), unpublished. 

4. D. Martello et al. 2003, 28th ICRC, Tsukuba. 

5. D. Heck et al., Report FZKA 6019 Forschungszentrum Karlsruhe. 

6. G. Di Sciascio and E. Rossi, Identification of showers with the core outside 
the ARGO- YBJ detector, these proceedings. 



